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ABSTRACT 



We have obtained 1.0 — 2.5^m spectra at R~600 of 14 disk dwarfs with spectral 
types M6 to L7. For four of the dwarfs we have also obtained infrared spectra at 
R~3000 in several narrow intervals. In addition, we present new L' photometry for 
four of the dwarfs in the sample, which allows improved determinations of their 
bolometric luminosities. While obtaining the photometry we resolved the L-dwarf 
Denis-P J 0205-1159 into an identical pair of objects separated by 0'.'35. The spectra, 
together with the published energy distribution for one other L5 dwarf, are compared 
to synthetic spectra generated by upgraded model atmospheres. Good matches are 
found for 2200>Tofj K>1900 (spectral types around M9 to L3), but discrepancies exist 
at Teff >2300 K (M8) and for Tog <1800 K (L4-L7). At the higher temperatures the 
mismatches are due to incompleteness in the water vapor opacity linelist. At the lower 
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temperatures the disagreement is probably due to our treatment of dust: we assume a 
photospheric distribution in equilibrium with the gas phase and neglect any diffusion 
mechanisms. We derive effective temperatures for the sample from the comparison 
with synthetic spectra and also by comparing our observed total intrinsic luminosities 
to structural model calculations (which are mostly independent of the atmosphere but 
are dependent on the unknown masses and ages of the targets). The two derivations 
agree to about 200 K except for the faintest object in the sample where the discrepancy 
is larger. Agreement with other temperature determinations is also ± ~200 K, except 
for the L7 dwarf. 



1. Introduction 

The last few years have seen rapid advances in both the observational data on low-mass stars 
and brown dwarfs and in the theoretical models of these objects. During this period we have 
been obtaining infrared spectra of a sample of halo and disk dwarfs approaching and below the 
stellar /sub-stellar boundary. This paper extends our spectroscopic study of M dwarfs presented 
in |Leggett et al. 2000a (hereafter LOO) to lower effective temperatures {T^s < 2500 K). The data 



are compared to synthetic spectra generated from an upgraded version of Allard & Hauschildt's 
NextGen model atmospheres ( Hauschildt, Allard &: Baron 199^ . 



In §2 we present the target sample and in §3 we describe the instruments used to obtain the 
spectroscopic data. In §4 we present imaging results: new L' photometry, confirmation of binarity 
of Denis-P J 0205-1159, and color-color diagrams. In §5 we describe the spectra and calculate 
integrated fluxes and bolometric magnitudes. §6 describes the models and the comparison process 
and the results of the comparison of the data to the synthetic spectra are presented in §7. In §8 
we present derivations of effective temperature based on structural models. Our conclusions are 
given in §9. 



2. The Sample 



The fourteen dwarfs whose spectra are reported here were selected from various color or 
proper-motion surveys — Delfosse et al. 1997| (the DENIS-P survey, sometimes known as DBD 
objects); prwin, McMahon Sz Reid 1991| (the BRI objects); Kirkpatrick, Beichman & Skrutskie 
[iMtI (the 2MASS survey); |Ruiz, Leggett k Allard T997| (Kelu-1); p?inney 1993| (the TVLM 
objects) — or from follow-up studies of known high proper-motion objects — Kirkpatrick, Henry 
|fc Irwin 1997 ; Monet et al. 199^ — or from studies of low-mass companions to nearby stars — 



Becklin &: Zuckerman 1988| (GD 165B). We also include for analysis the previously published 



energy distribution for the L dwarf SDSS 0539 found by the Sloan Digital Sky Survey (|Fan et 
al. 20001 , iLeggett et al. 2000b| ). The targets are listed in Table 1. We give discovery name, or 
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LHS or LP number (Luyten 1979), and/or Gliese or Gliese/Jahreiss number ( pliese fc Jahreis^ 



1991), for each dwarf. An abbreviated RA/Dec is also given to aid identification. The infrared 
instrumentation used for each object is given, using the configuration names from Table 2 (see §3). 

The spectral types in Table 1 are based on optical spectra and are taken from Kirkpatrick et 



al. 1995, 1999a, 2000, Kirkpatrick, Henry & Irwin 1997 and Fan et al. 200C. These classifications 



are all on the same system; we note the significantly earlier spectral type assigned to Denis-P J 



0205-1159 according to a classification system by pVlartfn et al. 1999| . The differences between the 
two classification systems are discussed further in §5.1. 



The kinematic populations have been taken from [Leggett et al. 1998| , [Leggett 1992| , except 
for Kelu-1, Denis-P J 0205-1159, Denis-P J 1058-1548 and Denis-P J 1228-1547 for which we have 



calculated space velocities using parallaxes and proper motions from Dahn et al. 2000 and radial 
velocities from Basri et al. 199^ , 2000. The space velocities of these four objects are low (U 20 to 
30 km/s, V —13 to —20 km/s, W 1 to 5 km/s) and imply young disk population. The kinematic 
classification schemes for young disk (YD), young/old disk (Y/0) and old disk (OD) are described 



m 



Leggett 1992. 



Detections of lithium in Kelu-1 ( [Ruiz, Leggett &: Allard 1997| , Basri et al. 1998), Denis-P J 
1228-1547 ( Pasri et al. 1998| , [Tinney, Delfosse fc Forveille 1997] ) and LP 944-20 ( flinney 1998D 
imply that each of these objects is substellar and at an age of ~1 Gyr. A non-detection for 
Denis-P J 1058-1548 ( Basri et al. 19981 ) implies a substellar mass if its age is ~1 Gyr or a mass 
at the stellar /substellar transition of 0.075 Mq if it is several Gyr. The age of GD 165B also is 
uncertain; its mass is near the stellar /substellar boundary, but it is likely to be a brown dwarf 



([Kirkpatrick et al. 1999bD . 

Denis-P J 1228-1547 was recently resolved to be a double system by jMartfn, Brandner & 



Basri 1999, using the Hubble Space Telescope. The binarity of Denis-P J0205-1159 was announced 



by Koerner et al. 1999 on the basis of K-band observations at Keck; in §4.2 we confirm this 



using JHKL' imaging data from the United Kingdom 3.8 m Infrared Telescope (UKIRT). 2MASP 
J0345432-I-254023 is suspected to be a double-lined spectroscopic binary ( pleid et al. 1999[ ). 
Indirect evidence for further multiplicity in our sample is discussed in §4.3. 



3. Spectroscopic Observations and Data Reduction 

The infrared spectra were obtained at UKIRT and at the University of Hawaii 2.2 m telescope, 
using some of the same spectrographs as LOO. Table 2 gives the observing dates, grating and 
slit information, and resolutions. The near-infrared data for SDSS 0539 were obtained using the 
configuration designated IRl in Table 1, except that the slit was wider (1'.'2), yielding resolutions 
of 25 A at J and 50 A at H and K. The lower resolution data cover the entire JHK bandpasses, 
from 1.0/im to 2.5/im. The medium resolution data cover narrower spectral regions: 1.19 — 1.30^m, 
1.58 — 1.67;um, and 2.02 — 2.18/im. The data are available in electronic form on request. 



- 4 - 



Data reduction was carried out in the usual way using the Figaro software package. The 
effect of the terrestrial atmosphere was removed by dividing by the spectra of nearby early-type 
stars, after removing hydrogen lines seen in their spectra. The shape of each infrared spectrum 
was corrected for the known flux distribution of the early-type star. The low-resolution spectral 
segments were individually flux-calibrated using previously reported or recently obtained JHK 
photometry. Each segment was integrated over the appropriate filter profile (UKIRT JHK) and 
the observed flux from Vega was integrated over the same profile; Vega was assumed to be zero 
magnitude at all wavelengths, and the target flux was scaled to match the broadband photometry. 
The higher resolution data were scaled to match the calibrated lower resolution data. 

Optical spectra were obtained from the literature for all the objects in the sample; the data 
sources are listed in Table 1. These data were flux-calibrated using the Cousins I filter and 
available I photometry, except for SDSS 0539 which was flux-calibrated using Sloan Digital Sky 
Survey photometry ( |Fan et al. 200C| ). 



4. Imaging Results 

4.1. New Photometry 

New infrared photometry was obtained for four of the objects in the sample. In 1999 February 
and March L' data were obtained for Denis-P J 1058-1548 and Denis-P J 1228-1547AB using 
IRCAM3 (pixel scale 0.28 arcseconds/pixel) at UKIRT. In 1999 September data were obtained 
for Denis-P J 0205-1159 (JHK and L') and 2MASP J0345432+254023 (L' only), using the same 
camera, which had been modified to a smaller pixel scale of 0.08 arcseconds/pixel (hereafter 
IRCAM/TUFTI). The results are presented in Table 3. 

Standards with a range of color were observed; the results indicated that no color term was 
needed to convert the V data taken with IRCAM/TUFTI to the estabhshed UKIRT (IRCAM3) 
L' system. However the J, H and K photometric systems are different due to significant differences 
in the filter set. The results in Table 3 are given in the new photometric system which we refer to 
here as the "UKIRT~UFTI(MKO-NIR)" system as the UFTI camera, with the new MKO-NIR 
filter set, is being used to define this system. Both UKIRT cameras, UFTI and IRCAM/TUFTI, 
are now configured with the new Mauna Kea consortium JHKL' filter set. Transformations 



between the established UKIRT system and the new MKO-NIR system are presented in [Hawarden 



et al. 2000 



4.2. Denis-P J 0205-1159 : an Identical Pair of L-Dwarfs 



Denis-P J 0205-1159, which was observed in excellent seeing conditions on 1999 September 
19 (UT) using the IRCAM/TUFTI camera, was easily resolved as a double system separated by 
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0.35ib0.03 arcseconds at a position angle of 77±4 degrees. Figure 1 shows the K image of the 
system. The JHKL' colors of each component are identical to within the measurement error. 
Proper motion measurements by Koerner et al. 199£ show that the object is a true binary. 
Those authors obtain a separation of 0.51±0.02 arcseconds and a position angle in 1997 July of 
106±5 degrees and, in 1999 January, 72ibl0 degrees. We suggest that our values of 0.35±0.03 
arcseconds and 77±4 degrees, based on higher resolution data (0.08 arcseconds/pixel cf. 0.15 
arcseconds/pixel) , are more accurate. The 0'.'35 separation translates to a physical separation of 
6.3 AU. This system should be monitored to obtain orbital information. 



4.3. Color Color Diagrams 

Table 4 lists distance moduli and VIJHKL' colors for the sample. The V and I photometry 



were obtained from |Dahn et al. 2000| for Denis-P J 0205-1159AB, Denis-P J 1058-1548, Denis-P 
J 1228-1547AB, 2MASP J0345432+254023 and Kelu-1. The I magnitude for SDSS 0539 is 
synthesized from the flux-calibrated spectrum. J, H and K magnitudes for this object are from 



Leggett et al. 2000b. The L' magnitude for GD 165B is taken from [Jones et al. 1996 . Otherwise 



the photometry is from this work or taken from the compilations by [Leggett 1992 , Leggett et al. 



1998. The infrared photometry is on the established UKIRT (IRCAM3) system. We have averaged 



the JHK values for Denis-P J 0205-1159AB presented here with the earlier values in Leggett et 



al. 1998| after transforming the new JHK to the old photometric system using the known filter 



profiles. The two sets of values agree well. 

Figures 2-4 show various VIJKL' color-color diagrams; although similar diagrams are shown 



in Leggett et al. 1998 the ones shown here use improved or additional data at V, I and L', 
and are useful references as clear trends are seen. Object labels in these and following plots 
are abbreviated: Denis-P J0205-1159AB, Denis-P J 1058-1548, Denis-P J 1228-1547AB become 
DNS 0205, 1058 and 1228 respectively; 2MASP J0345432-h254023 is labelled as 2M 0345; TVLM 
513-46546 is shown as T 513-46546. 

Figure 3 shows a synthetic color-temperature sequence for solar-metallicity and log{g) =5.0, 
generated from Ames-Dusty-2000 models (see §6) by convolving the calculated energy distributions 
with the known UKIRT filter profiles. There are some discrepancies between the calculated 
and observed colors due to known inadequacies in the models which we discuss in §6 and §7. 
Nevertheless this comparison provides the temperature range required for the spectral analysis 
model grid (approximately 1800 K to 2300 K). 



Figure 5 shows Mj:J— K with isochrones from phabrier et al. 2000[ for ages 5 Gyr and 0.1 Gyr. 
These isochrones use an earlier very similar Dusty version of the NextGen model atmospheres. 
The agreement is poor for the hotter dwarfs with Mj <12, and not good for the coolest objects 
Denis-P J 0205-1159 and Denis-P J 1228-1547. The discrepancies reflect known inadequacies in 
the atmospheric models used to generate the colors (§6) rather than problems with the structural 
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models. 

As mentioned in §3, LP 944-20, Kelu-1 and Denis-P J 1228-1547AB are young substellar 
objects, and GD 165B (L4) and Denis-P J 1058-1548 (L3) are substellar or stellar /substellar 
transition mass objects, depending on their age. Kelu-1 is superluminous but no companion is 
seen in HST images ( [Martm, Brandner fc Basri 1999| ). Martin et al. 1999| point out that if Kelu-1 



is a single object its luminosity implies an age of <0.1 Gyr. [Basri et al. 2000| find that the object 
is a very rapid rotator (vsini = 60km/s) and so it may indeed be very young, if the correlation 
seen for hotter low-mass stars between youth and rapid rotation holds here. Nevertheless it is 
important to monitor it for radial velocity variations. 2MASP J0345432-I-254023, which also 
appears superluminous, is suspected to be a spectroscopic binary (see §2). Both Denis-P J 
0205-1159 and Denis-P J 1228-1547 are known to consist of pairs of very similar objects. The 
high frequency of binaries is not surprising in an effectively magnitude-limited sample such as 
this. Note that the evolutionary models imply that each component of the ~L5 dwarfs Denis-P J 
0205-1159 and Denis-P J 1228-1547 must be substellar even for an age of several Gyr. 



5. Spectroscopic Results 

5.1. Spectroscopic Sequences 

Figure 6 shows a representative set of R~600 spectra covering the temperature range of our 
sample, which corresponds to spectral types dM7 to around dL7. The obvious features to note 
are: the CO bands at 2.3-2.4/xm; the water bands at 1.4^m, 1.8/xm and longward of 2.4/im; the 
FeH bands at 0.99^m and 1.2//m; and the KI doublets near l.lS^m and 1.24/um. Note that the 
TiO bands, which are used to classify the M dwarfs, are not apparent or very weak in the L dwarf 



spectra. A pseudo-continuum approach has been used by Martin et al. 1991; to classify the L 
dwarfs from optical spectra. They find a good correlation between this scheme and the strength 
of the water band at 1.4^m. The trend of strengthening water absorption with later spectral type 



is seen in Figure 6. Kirkpatrick et al. 1999a use ratios of various optical absorption features to 
classify the L dwarfs. Both schemes roughly agree up to L5 at which point the Martin classes 
become earlier than Kirkpatrick's (so that, for example, Denis-P J 0205-1159 is classified L5 in the 
Martin scheme and L7 on Kirkpatrick's). 

We note that the shape and depth of the water bands are calculated to be strongly affected 
by the presence of dust, which affects the emergent spectra for T^g < 2500 K, or spectral types 
dM7 and later. With no dust in the atmosphere, the bands are predicted to become deeper with 
increasingly steep wings as temperature decreases. When dust is present the atmosphere is heated 
in the line forming region and the bands become shallower and wider. It is likely that as infrared 
spectra are obtained for a larger sample of L dwarfs, more scatter will be seen in the relationship 
between the depth of the water bands and effective temperature or (optically derived) spectral 
type, as dust properties will vary from object to object (due to e.g. rotation, metallicity or age 
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difFerences) . 



5.2. Integrated Fluxes and Bolometric Corrections 

Table 4 gives integrated fluxes for the sample, expressed as flux at the Earth, apparent 
bolometric magnitude and intrinsic stellar luminosity. They were obtained by integrating the 
infrared spectra over wavelength and adding estimated flux contributions at shorter and longer 
wavelengths. Each shorter wavelength contribution was approximated as a linear extrapolation 
to zero flux at zero wavelength. The contributions at wavelengths beyond 2.4/im were calculated 
by deriving the fluxes at L' using an effective wavelength approach, interpolating the spectrum 
from the end of the K-band spectrum to the L' wavelength, and assuming a Rayleigh-Jeans tail 
beyond L' (model calculations imply that the Rayleigh-Jeans assumption is good to ~1% at 
these temperatures). Although it is now known that methane absorption at 3.3/im is seen in late 
L-dwarfs ( Noll et al. 2000| ), the feature is weak and will not significantly effect the determination 



of total flux for this sample. For SDSS 0539 we adopted an L' magnitude based on J— K, the 
K magnitude, and Figure 4. We note that for dwarfs with 2700 K > Tcs > 2000 K (spectral 
types later than dM5.5) the flux beyond 2Anm makes up >20% of the total flux and hence L' 
measurements are crucial for an accurate determination of luminosity. The uncertainties are 5%, 
0.05 mag and 0.02 dex in total flux, bolometric correction and logio-^^/i©, respectively. 

Figure 7 plots the K-band bolometric correction against I— K and Figure 8 shows JiCx 
against J— K. We have included the results from Leggett et al. 199^ and LOO. For the non-halo 



(disk and likely disk) dwarfs only, the relationship between K-band bolometric correction and 
color is represented by the cubic polynomials: 

BCk = 1.307 + 0.9147(1 - K) - 0.1592(1 - Kf + 0.01054(1 - Kf 

BCk = -0.31 + 5.124(J -K)- 2.031(J - Kf + 0.13877(J - Kf 

for 1.9 < I— K < 5.5 and 0.75 < J— K < 1.60. These relationships are indicated by the solid lines 
in Figures 7 and 8, and supercede that given in LOO. 



6. Models and Synthetic Spectra, and Comparison Process 

The Ames-Dusty-2000 models used for this work were calculated as described in LOO, and 
are based on the Ames H2O and TiO line lists by Partridge and Schwenke 1997 and Schwenke 
We stress that large uncertainties persist in these opacities for the temperature range of 



1998 



this work (see Allard, Hauschildt &: Schwenke 2000|) . The models have been upgraded with (i) 
the replacement of the JOLA (Just Overlapping Line Approximation) opacities for FeH, VO and 
CrH by opacity sampling line data from Phillips et al. 1987 and R. Freedman (NASA-Ames, 
private communication), and (ii) the extension of our database of dust grain opacities from 10 to 
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40 species. These models and their comparisons to earher versions are the subject of a separate 
pubhcation (Allard et al., in preparation). 

The fitting of the synthetic to the observed spectra was done using an automatic IDL program. 
First, the resolution of the synthetic spectra was degraded to that of each individual observed 
spectrum and the spectra were normalized to unit area for scaling. The comparison was done 
using a model atmosphere grid that covers the range 1500 K< Tcs < 3000 K, 3.5 <log{g) < 6.0 
and [M/H]= —1.0,-0.5 and 0.0, with a total of 162 model atmospheres. For each observed 
spectrum we then calculated a quality function q, similar to a value, for the comparison with all 
synthetic spectra in the grid. In order to avoid known problematic intervals in either the observed 
or the synthetic spectra, the wavelength range 1.5-1.7//m was excluded from the comparison. Two 
particularly diagnostic wavelength ranges, 0.7-1.4//m and 2.0-2.5/im, were given 5 times higher 
weights than the remaining spectral ranges. We selected the models that resulted in lowest 3-5 
q values as the most probable parameters range for each individual star. The "best" value was 
then chosen by visual inspection. This procedure allows a rough estimate of the uncertainty in the 
stellar parameters. Note that it does not eliminate systematic errors in the stellar parameters due 
to missing or incomplete opacity sources. 

This set of models relies on the assumption that the dust is in equilibrium with the gas phase 
at any depth in the atmosphere. This corresponds to a limiting case, discussed also by Chabrier 
let al. 2d00| in the context of brown dwarf evolution, which could be met if the mixing of material 
(by convection, rotation or other processes) prevents the diffusion of the dust. Of course this 
situation may not be met at all layers of the atmosphere and/or at all effective temperatures 
since, for example, processes such as convective mixing retreat from the upper atmosphere as Teg 
decreases. The current analysis will help determine the parameter range of stars over which this 
approximation is valid. 



7. Results of Comparison of Data and Atmospheric Models 

7.1. Low Resolution Spectra; Derived Temperatures and Diameters 

Figure 9 shows examples of our best model fits to the observed spectral energy distributions. 
Table 5 lists the best-fit model parameters, Tgfj, log(g) and [m/H], for each dwarf, sorted by 
spectral type (on the Kirkpatrick scheme). We also give diameters for the fourteen dwarfs with 
parallax measurements, derived in two ways — one from the scaling factor required to match the 
calculated surface flux to that observed at the Earth, and the other from the measured luminosity 
and derived effective temperature, via the Stefan-Boltzmann law. We return to the diameter 
measurements in §8. 

The synthetic spectra fit the energy distributions for early-L dwarfs with Tgfj ~1900 K quite 
well, but there are obvious discrepancies for the hotter and cooler dwarfs. Discrepancies for late-M 
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dwarfs are also apparent in LOO — the predicted water bands are too deep and there is too much 
flux predicted at the peaks of the J and H bands (in other words, the water vapor opacity profile 
is too cold). We attribute this to remaining problems in the calculated water opacity database, 



which is known to be incorrect for higher temperature transitions (Allard, Hauschildt & Schwenke 



2000). There are also discrepancies seen in the fits to the two latest spectral types in our sample, 
Denis-P J 1228-1547AB and Denis-P J 0205-1159AB. Although the shapes of the water bands 
are better matched for these objects, the calculated fluxes are too high at the peaks of the J and 
H bands and the models do not reproduce the shape of the K-band peaks. This indicates an 
inaccuracy in the atmospheric structure which appears to be too hot in the deeper layers from 
which the J- and H-band flux emerges. These discrepancies reveal the increasing inadequacy 
of the dust equilibrium assumption (or neglect of grain diffusion mechanisms) in the model 
photosphere construction as T^s decreases. While models which incorporate gravitational settling 
by diffusion are not yet available, our result demonstrates the importance of such mechanisms in 
the L dwarf temperature range. Our determination that Denis-P J 1228-1547AB is metal-poor 
may only reflect problems in the Dusty models at this Tgfj — it may in fact have solar metallicity 
(see §7.2). Note that despite the improved treatment of FeH, Figure 9 shows that better line data 
are still needed for this important opacity source. Note also that investigations of the effect of 



using different water opacities are presented in Allard, Hauschildt fc Schwenke 200C. 



For this particular model grid comparison the range of T^q allowed by the spectral energy 
distribution fitting is constrained to ±75-100 K, however this does not include systematic effects. 
There is further discussion of the errors in our Tefj values in §8 and §9. Surface gravity also affects 
the overall shape of the energy distribution. The best-fit models have log(g') of 5.5 or 6.0, except 
for the apparently metal-poor and slightly lower gravity Denis-P J 1228-1547AB. Given the 
remaining uncertainties tied to the water opacity and treatment of dust, we estimate the accuracy 
of our determinations of log{g) to be 0.5 dex. The structural models of phabrier et al. 2000| yield 
log(g) ?a5.3 for objects in the observed temperature range with ages 1 — 5 Gyr. 



7.2. Medium Resolution Spectra 

Figure 10 shows the spectra we obtained at R~3000. Only the J- and K-band regions are 
shown; the spectra of a portion of the H-band contain no obvious features. Also shown are the 
best-fit model synthetic spectra, based on the match to the lower resolution data but shown at 
higher resolution. We have rescaled the synthetic spectra to match the local pseudo-continuum. 

The dominant narrow features in the J-band are the potassium doublets. The model spectra 
generally reproduce the shape and depth of the doublet, as well as the structure to the blue of 
the doublet. The pronounced asymmetry in the K I doublet for Denis-P J 0205-1159AB is real, 
but is not understood at this time. Most of the fine structure in the K-band is real and is due to 
water absorption; the synthetic spectra reproduce the observed structure well even at this detailed 
level. Experiments with varying log(g) and [m/H], while keeping the effective temperature fixed, 
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showed that at R~3000 the hne profiles did not constrain gravity or metahicity significantly better 
than achieved with the R~600 data (but they did support the parameters derived from the other 
data). For the Denis L dwarfs the J-band data were more diagnostic than the K-band, but only 
constrained log{g) to values between 5.0 and 5.5, and for Denis-P J 1228-1547AB showed that a 
metallicity of [m/H]= was slightly preferable to the [m/H]= —0.5 implied by the low resolution 
fit. 



8. Comparison of Data with Structural Models 

Figure 11 shows the effective temperature implied by our atmospheric model fitting versus 
diameter derived by scaling the model flux to that observed at the Earth. We also show calculated 



sequences for ages 0.1 Gyr and 5 Gyr from Chabrier et al. 2000. There are discrepancies between 



our data points and the structural model at both the hot and cool ends of our sample. The 
diameters derived here are subject to errors in the parallax measurements as well as in the 
calculated surface flux. The variation in flux levels between the next-best-fit model and the 
best-fit model imply a typical total error in the diameter of 10 — 15%. The error in the diameter 
determined from luminosity and Tefj is dominated by the uncertainty in Tefj and is typically 15% 
if the error in T^s is ~100 K. These are the errors shown by the error bars in Figure 11. 

However as discussed above there are large uncertainties remaining in the opacity treatment 
for the temperature range of this work (see Allard, Hauschildt fc Schwenke 2000| ), and in the 



treatment of dust — the assumption here that the dust is in equilibrium with the gas phase at any 
depth in the atmosphere corresponds to a limiting case, discussed also by phabrier et al. 200C . It 



is possible that systematic errors exist in our derivations of effective temperature and diameter, 
despite our exclusion of the most suspect wavelength regions in the model fitting. 

The theoretical determinations of radii and luminosity for low-mass stars and brown dwarfs 
are thought to be reasonably robust and reasonably independent of the surface atmosphere used 



for the structural model. This is demonstrated by the fact that the radii calculated by Burrows et 



al. 19971 for brown dwarfs of mass 0.03—0.07 M0 with 1300 < T^s K< 2600, and age 0.3—3.0 Gyr, 



range (only) from 5.8e7 m to 7.6e7 m, in close agreement with those derived by Chabrier et alj 



2000| of 6.0 — 7.6e7 m. Figure 12 plots Teg versus logioL/L© for objects ranging in mass from 



0.03 Mq to 0.10 Mq and aged 0.1 Gyr to 10 Gyr. These are a reasonable to generous range in 



parameters for our observed range in luminosity. The sequences shown are from Chabrier et al 



2000, but one mass sequence from Burrows et al. 1997 is also shown to demonstrate the good 



agreement between evolutionary models. We have used this figure to derive effective temperatures 
for the targets in the sample with known distance (and hence luminosity) and these values are 
listed in the last column of Table 5. We assume all objects are single except for Denis-P J 
1228-1547AB and Denis-P J 0205-1159AB, for which we assume two identical components. Note 
that the older kinematic population class for LHS 429 would suggest that the effective temperature 
lies at the high end of the tabulated range for this object. 
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There is an indication that our effective temperatures derived from fitting the spectral energy 
distributions with the synthetic spectra may be underestimated for the hotter stars in our sample, 
with spectral type around M7, and overestimated for the cooler mid- to late-L dwarfs. A hotter 
scale for the late-M dwarfs would agree better with that derived in Leggett et al. 2000b| (i.e. 
Tefj ~ 2600 K for type dM6). The agreement between these two approaches — the spectral 
synthesis and luminosity comparison — is ~200 K, except for the coolest object in the sample, 
Denis-P J 0205-1159AB, where the discrepancy is ~400 K. Better determinations of effective 
temperature for objects this late will have to await models with improved treatment of dust. 



9. Discussion of the L Dwarf Temperature Scale and Conclusions 

Our derived temperatures agree with other derivations of T^q for early-L type dwarfs within 
the likely errors. Pavlenko et al. 200C| use both the NextGenDusty and Tsuji 200C| atmospheric 



structures with their own spectral synthesis code to fit red spectra; they have to include an 
unknown additional opacity source, in the form of a power-law, to match the overall spectral 
shape because they do not include dust in solving for chemical equilibrium, or account for the 
opacity of known dust species (note that this is inconsistent with their adopted atmospheric 
structure). They derive T^s for four dwarfs in common with us: for BRI 0021-0214 and Kelu-1 
their values are higher than our spectral synthesis values by 100 K; for Denis-P J 1228-1547AB 
and Denis-P J 0205-1159AB their values are 200 K and 700 K lower. Their value of 1200 K for 
Denis-P J 0205- 11 59 AB does not seem likely, as we would then expect to see methane features in 
the and K-band regions, which are not seen. 
I I 

Basri et al. 2000 use Ames-Cond NextGen model fits to rubidium and cesium lines and 
derive effective temperatures within 100 K of our spectral synthesis values for Kelu-1, Denis-P 
J 1058-1548 and Denis-P J 1228-1547AB; they derive a value 150 K lower than ours for the 
coolest object in common, Denis-P J 0205-1 159AB, and a value 300 K higher for the M9.5 dwarf 
LP 944-20. The Ames-Cond models include the formation of dust but neglect its opacity (i.e. 
assume the dust has settled below the photosphere), which corresponds to a limiting case which 
complements that assumed here of the dust in equilibrium with the gas phase. It is expected that 
dust clouds will sink lower in the atmosphere with decreasing effective temperature (in the absence 
of replenishment of material in the upper photosphere), hence the Cond models are likely to be 
more correct than the Dusty models for the latest L-dwarfs. The Basri et al. lower temperature 
of 1750 K for Denis-P J 0205-1159AB is likely to be more accurate than that derived here, but 
the opposite is true at hotter temperatures and so our lower Teg for LP 944-20 is likely to be more 
accurate than their value. Given the discrepancies seen in the model fits to the observed spectral 
energy distributions for the coolest objects in our sample, it would appear that gravitational 
settling of dust becomes important at Tefj < 1800 K. 



Based on species condensation arguments, Kirkpatrick et al. 2000 argue for an effective 
temperature around 2000 K for LO spectral types, which is in agreement with this work. They 
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also advocate (based on luminosity) T^s ~1300 K for the L8 dwarf (by the Kirkpatrick scheme) 
GL 584C. However the discovery of three dwarfs optically redder than L8 and showing both CO 
and CH4 in their spectra ( Leggett et al. 2000b| ) suggests that these new later-type objects have 



Teff -^1300 K and by implication L8 spectral types have T^s >1300 K. A hotter temperature scale 
also is supported by the preliminary model fits to L-band methane features seen in mid-L dwarfs 
reported by Noll et al. 2000| . However all of these results can be estimates only as they either rely 



on the correlation between the appearance or disappearance of various species, which depends on 
the temperature of the line-forming region and not on the effective temperature, or they rely on 
what are acknowledged to be very preliminary atmospheric models in this temperature range. 

In conclusion, the best current estimates are that the M dwarfs cover a range in T^s of 
3700—2100 K, the L dwarfs 2000~~1500 K and the known methane (T) dwarfs ~1300— -800 K. 
A definitive temperature scale for the L dwarfs will not be possible until model atmospheres which 
include gravitational settling are available; such work is in progress by the modellers of this group 
and others. We note that any object with a spectral type L5 or later must be substellar based 
on evolutionary models, and late-M to early-L dwarfs may be substellar depending on their age. 
The discovery of binary L-dwarf systems offers a chance to further improve our knowledge of the 
fundamental parameters for this recently discovered low-mass population of our Galaxy. 
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Fig. 1.— K band image of Denis-P J 0205-1159 
to the left. 



The pair are separated by 0'.'35; North is up, East 
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Fig. 2. — V— 1:1— K diagram; filled symbols are this work, open symbols are dwarfs from Leggett 
et al. 1998 , Leggett 1992| (on the UKIRT JHK system). Symbol shapes represent kinematic 
populations: triangles — disk, circles — unknown. Typical uncertainties are shown at lower right. 
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Fig. 3. — I— J:J— K diagram using the same symbols as in Figure 2. Spectral types based on I— J 
are shown for guidance only — for the L dwarfs there is significant variation of color with type. The 
solid line is the Dusty sequence for [m/H]=0 and log{g) =5.0 as a function of effective temperature. 
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Fig. 4. — J— K:K— L' diagram, using the same symbols as in Figure 2. 
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Fig. 5. — Mj:J— K with symbols as in Figure 2. The location of the identical components of Denis-P 
J 0205-1159 and Denis-P J 1228-1547, if they were resolved, are indicated by the dashed line and 
lower triangles. The isochrones for 0.1 Gyr and 5.0 Gyr are from Chabrier et al. 2000. 
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Fig. 6. — Spectral sequence for late-M to L dwarfs where the spectra have been normahsed to the 
flux at 1.2fim and offset. Horizontal dotted lines indicate zero flux levels. Denis-P J 0205-1159 is 
classifled L5 in the Martin scheme. 
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Fig. 7. — Bolometric correction at K versus I— K. Filled symbols are this work, open are from 
Leggett et al. 1996, 2000a. Symbol shapes represent kinematic populations: squares — halo, 
triangles — disk, circles — unknown. The solid line is the empirical fit to the disk objects. 
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Fig. 8. — Bolometric correction at K versus J— K. Symbols are as in Figure 7. The solid line is the 
empirical fit to the disk objects (circles and triangles). 
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Fig. 9. — Examples of best synthetic spectra fits to late-M and L dwarfs, where the models are 
shown as dotted lines and Tcff/log((7)/[m/H] are given. The spectra have been normalised to the 
flux at 1.2/im and offset. Horizontal dotted lines indicate zero flux levels. Denis-P J 0205-1159 is 
classified L5 in the Martin scheme. 
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Fig. 10. — Higher resolution data with synthetic spectra for the best fit model as determined from 
the low resolution spectral energy distributions. The spectra have been normalised and offset for 
clarity. Much of the structure in the K-band is real and due to water. 
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Fig. 11. — Diameter (derived by scaling) as a function of T^s, from the model atmosphere 
comparison. Dashed lines are structural model calculations from Chabrier et al. 2000| for 0.1 Gyr 
and 5.0 Gyr. The location of the identical components of Denis-P J 0205-1159 and Denis-P J 
1228-1547, if they were resolved, are indicated by the dashed lines to smaller diameters. 
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Fig. 12. — Luminosity as a function of Tefj for constant mass, from the structural models of Chabrier 
et al. 2000| . Masses are indicated in units of M0 and evolution proceeds right to left; dotted lines 



indicate 0.1, 1.0 and 10.0 Gyr. The circles indicate the 0.04 Mq sequence from Burrows et al. 1997 
for comparison. 
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Table 1. Dwarfs in the Sample 



LHS 


Other 


RA/Dec 


Spectral*^ 


Kinematic 


Instrumentation' 


Number 


Name 


1950 


Type 


Population 


Infrared 


Optical 




BRI 0021-0214 


0021-02 


M9.5 


Y/0 


ixvl 


1 
1 




Denis-P J 0205-1159AB 


0203-12 


L7 


YD 




1 




LP 944-20 


0337-35 


M9 


YD 


IRl 


2 




2MASP J0345432+254023 


0342+25 


LO 




IRl 


3 




SDSS 0539 


0539-00 


L5 




4 


5 


zUdo 


LP DDD-9 


0851—03 


M9 


YD 


IR1,IR0 


6 




Denis-P J 1058-1548 


1058-15 


L3 


YD 


IR1,IR0 


1 




Denis-P J 1228-1547AB 


1228-15 


L5 


YD 


IR1,IR0 


1 




Kelu-1 


1305-25 


L2 


YD 


IR3 


7 




GD 165B 


1422-F09 


L4 




IR3 


6 


2930 


LP 98-79 


1429+59 


M6.5 


Y/0 


IR4 


8 


3003 


LP 914-54 


1453-27 


M7 


Y/0 


IR4 


9 




TVLM 513-46546 


1458+23 


M8.5 


Y/0 


IR3 


9 


429 


Gl 644C 


1652-08 


M7 


OD 


IR4 


9 


3406 


LP 229-30 


1841+40 


M5.5 


Y/0 


IR4 


10 



^Spectral types quoted on the Kirkpatrick scheme, note that Denis-P J 0205-1 159 AB is classified 
L5 on the Martin scheme, see text 

^Published spectra taken from: 

1 Tinney et al. 1998 resolution 7A 

2 Tinney & Reid 1998 resolution 0.4A 

3 Kirkpatrick, Beichman & Skrutskie 1997 resolution 12A 

4 Leggett et al. 2000b resolution 25— 50A 

5 Fan et al. 2000 resolution 14A 

6 Kirkpatrick et al. 1995 resolution 18A 

7 Ruiz et al. 1997 resolution 17A 

8 Kirkpatrick 1992 resolution 18A 

9 Kirkpatrick et al. 1995 resolution 12A 

10 Reid, Hawley & Gizis 1995 resolution 8 A 
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Table 2. Instrumentation for New Spectroscopic Data 



Configuration Instrument Telescope Date Grating Slit Resolution A 

Name YYYYMM lines/mm Width" Z J H K 



IRO 
IRl 



IR3 
IR4 



CGS4 
CGS4 



CGS4 
KSPEC 



UKIRT 
UKIRT 



UKIRT 
UH 88" 



199801 
199904 
199801 
199711 
199704 
199407 



150 
40 



75 
60 



1.2 
0.6 



1.2 
1.2 



4 6 6 
17 35 35 



•• 16 30 30 
20 25 30 40 



Table 3. New Photometry 



Name 


J (error) H (error) K (error) 


L' (error) 




UKIRT-UFTI(MKO-NIR) 


UKIRT 



Denis-P J 0205-1 159AB 
2MASP J0345432+254023 
Denis-P J 1058-1548 
Denis-P J 1228-1547AB 



14.43 (0.04) 13.61 (0.03) 12.99 (0.03) 



11.44 (0.10) 
12.01 (0.10) 
11.62 (0.10) 
11.42 (0.10) 
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Table 4. Colors (on the Cousins and UKIRT-IRCAM3 systems) and Fluxes 



Name 


M 




V- 


-I 


I 


I- 


-J 


J- 


H 


H- 


K 


K 


K-: 


L' 


Flux 


b 




BC 


d 

K 


log 10" 


































W/m 


2 












BRI 0021-0214 


-0, 


.42^ 






15, 


,10 


3, 


,30 


0, 


,76 


0, 


,46 


10.57 


0, 


,79 


8.17c- 


14 


13, 


.74 


3, 


.17 


-3.43 


Denis-P J 0205-1 159 AB 


-1. 


.28^ 






18, 


.44 


3, 


.94 


0, 


,98 


0, 


.51 


13.01 


1, 


,57 


8.60e- 


15 


16, 


.19 


3, 


,18 


-4.06 


LP 944-20 


1. 


.511 






14, 


.16 


3, 


.48 


0, 


,70 


0, 


.45 


9.53 


0, 


,81 


2.20e- 


13 


12, 


.67 


3. 


,14 


-3.77 


2MASP J0345432+254023 


-2, 


.182 


4, 


.66 


17, 


,36 


3, 


,45 


0, 


,77 


0, 


,44 


12.70 


0, 


,69 


1.12e- 


14 


15 


.90 


3, 


,20 


-3.59 


SDSS 0539 










17, 


,67 


3, 


.73 


0, 


.97 


0, 


.53 


12.44 






1.34e- 


14 


15, 


.70 


3, 


,26 




LHS 2065 


0, 


.35* 


4, 


.36 


14, 


.44 


3, 


.22 


0, 


,80 


0, 


.48 


9.94 


0, 


,75 


1.45e- 


13 


13 


.12 


3. 


,18 


-3.49 


Dcnis-P J 1058-1548 


-1, 


.22^ 






17, 


,80 


3, 


,61 


0, 


,99 


0, 


,59 


12.61 


0, 


,99 


l.lOe- 


14 


15, 


.92 


3, 


,31 


-3.98 


Denis-P J 1228-1547AB 


-1, 


.292 






18, 


.21 


3, 


.83 


1, 


,05 


0, 


.57 


12.76 


1, 


,34 


9.94e- 


15 


16 


.03 


3, 


,27 


-4.00 


Kelu-1 


-1, 


.442 


4, 


.84 


16, 


.94 


3, 


.65 


0, 


,92 


0, 


.55 


11.82 


1, 


,04 


2.37e- 


14 


15, 


.09 


3. 


,27 


-3.56 


GD 165B 


-2, 


.49* 






19, 


,16 


3, 


,45 


1, 


,03 


0, 


.55 


14.13 


1, 


,20 


2.82c- 


15 


17 


.40 


3, 


,27 


-4.06 


LHS 2930 


0, 


.08* 


4, 


.57 


13, 


.31 


2, 


.54 


0, 


,65 


0, 


.38 


9.74 


0, 


,52 


2.05e- 


13 


12, 


.74 


3, 


,00 


-3.23 


LHS 3003 


0, 


.97* 


4, 


.52 


12, 


.53 


2, 


.51 


0, 


,65 


0, 


.42 


8.95 


0, 


,52 


4.14e- 


13 


11, 


.98 


3, 


,03 


-3.28 


TVLM 513-46546 


0, 


.04^ 






15, 


,09 


3, 


,29 


0, 


,68 


0, 


.40 


10.72 


0, 


,69 


7.43e- 


14 


13 


.85 


3, 


,13 


-3.65 


LHS 429 


0, 


.95* 


4, 


.56 


12, 


.24 


2, 


.39 


0, 


,63 


0, 


.38 


8.84 


0, 


,47 


4.69e- 


13 


11, 


.84 


3, 


,00 


-3.22 


LHS 3406 


-0. 


.75* 


4. 


.36 


13, 


.87 


2, 


.56 


0, 


,61 


0, 


.35 


10.35 


0, 


,57 


1.24e- 


13 


13, 


.29 


2. 


,94 


-3.12 



^Distance moduli from: 

1: Tinney 1996 

2: Dahn et al. 2000 

3: Yale General Gatalogue (van Altena et al. 1995) 
''Flux at the Earth 

'^Apparent bolometric magnitude; adopting L© = 3.86e26 W and MioIq = 4.75 then: mboi = —2.5 x logio{flux) — 18.978 
"^BCk = ruboi- K 

''Intrinsic luminosity, adopting L© = 3.86e26 W then: IoqioL/Lq = logio{flux) — 2 x logioT^ + 7.491 
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Table 5. Derived Parameters for the Sample 



Spectral Name By Spectral Synthesis By Structural Model 

Type Teff K/log(£?)/[m/H] Diameter 10^ m Teff K 

by Scaling by L,Tef[ 



M5.5 


LHS 3406 


2300/ 


6.0/ 


0.0 


2.38 


2.42 


2550-2700 


M6.5 


LHS 2930 


2300/ 


6.0/ 


0.0 


2.09 


2.14 


2400-2650 


M7 


LHS 3003 


2200/ 


6.0/ 


0.0 


2.20 


2.20 


2400-2650 


M7 


LHS 429 


2300/ 


6.0/ 


0.0 


2.13 


2.16 


2450-2650 


M8.5 


TVLM 513-46546 


2200/ 


6.0/ 


0.0 


1.45 


1.44 


2000-2250 


M9 


LP 944-20 


2100/ 


6.0/ 


0.0 


1.43 


1.38 


1850-2100 


M9 


LHS 2065 


2100/ 


6.0/ 


0.0 


1.97 


1.90 


2150-2400 


M9.5 


BRI 0021-0214 


2100/ 


6.0/ 


0.0 


2.12 


2.03 


2200-2450 


LO 


2MASP J0345432+254023(SB?) 


2000/ 


6.0/ 


0.0 


1.92 


1.87 


2050-2350 


L2 


Kelu-1 


1900/ 


5.5/ 


0.0 


2.10 


2.14 


2100-2350 


L3 


Denis-P J 1058-1548 


1900/ 


5.5/ 


0.0 


1.28 


1.32 


1700-1950 


L4 


GD 165B 


1900/ 


5.5/ 


0.0 


1.14 


1.20 


1650-1850 


L5 


SDSS 0539 


1900/ 


5.5/ 


0.0 








L5 


Denis-P J 1228-1547AB 


1800/ 


5.0/- 


-0.5 


1.42 


1.44 


1450-1650^ 


L7 


Denis-P J 0205-1159AB 


1900/ 


5.5/ 


0.0 


1.11 


1.20 


1400-1600' 



^Assuming equal contribution to luminosity from each component of binary 



